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Abstract. Terrestrial small mammal community composition was examined in the Volcanoes National Park (VNP), Rwan- 
da with respect to variation in habitat type and altitude. Trapping resulted in the capture of 220 individuals including 
eight species of rodents, three species of shrews and one mongoose. Of the species captured, Praomys degraaffi is vul- 
nerable and Sylvisorex vulcanorum 1s near threatened (IUCN 2012). Six species (Hylomyscus vulcanorum, Mus bufo, 
Praomys degraaffi, Sylvisorex vulcanorum, Lophuromys woosnami and Tachyoryctes ruandae) are endemic to the Alber- 
tine Rift. Species richness and diversity varied significantly among the different habitat types as they increased with el- 
evation up to the middle altitudes (2860—3255m) and then declined with increasing elevation. Altitude accounted for 
20 % of the variation in species diversity. Species variation in habitats and altitude was related to environmental factors. 
Endemic species were found mainly in low and middle attitude habitats; thus these habitat types are important for con- 
servation of small mammals at Volcanoes National Park. 
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INTRODUCTION 


Baseline data are essential for biodiversity monitoring, es- 
pecially during this era of anthropogenic and climatic 
change. However, the lack of accurate data on the status 
of free-ranging wildlife populations may limit the effica- 
cy of monitoring programs. The African montane primary 
forests of the Albertine Rift are home to many endemic 
mammal species (Hutterer et al. 1987; Kerbis Peterhans 
et al. 1998; Kaleme et al. 2007). The Volcanoes Nation- 
al Park, gazetted as Parc National des Volcans in 1925 
(hereafter VNP) is an afromontane forest rich in mammal 
species and endemism and is continuous with adjacent vol- 
canoes in the Democratic Republic of Congo and Ugan- 
da. Many small mammals exhibit narrow habitat prefer- 
ences and the distribution and abundance of small mam- 
mals may be indicators of ecosystem health (Golley et al. 
1975). Habitats with low vegetation diversity house rel- 
atively few species (Kaleme et al. 2007), and altitudinal 
variation also influences small mammal diversity (Kerbis 
Peterhans et al. 1998). 

The VNP is thought to support more than 25 small mam- 
mal species (Gyldenstolpe 1928; Hutterer et al. 1987). 
However, few studies on the ecology and community char- 
acteristics of small mammals (rodents and shrews) have 
been done in the VNP region (Gyldenstolpe 1928; Hut- 
terer et al. 1987). While a considerable amount of work 
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has been carried out on the vegetation, birds and large 
mammals of VNP, including the critically endangered Go- 
rilla beringei beringei Matschie, 1914 (Plumptre 1991; 
Fischer & Hinkel 1992; Robbins et al. 2001; Owiunji et 
al. 2005), small mammals remain poorly documented. The 
only scientific publication on small mammals in the VNP 
was written over 80 years ago (Gyldenstolpe 1928) when 
20 species were noted (data included in Table 1), but de- 
tailed small mammal diversity and distribution patterns 
were not determined. Since then, there has only been one 
unpublished study by Kajonjoli in 1993, rendered incom- 
plete by 1990's Rwandan civil war and the unpublished 
report (data included in Table 1) on the Ugandan slope of 
Mgahinga Gorilla National Park (MGNP) by Kerbis Pe- 
terhans & Austin (1996). 

Information on small mammal distribution and abun- 
dance across diverse habitats and along an elevational gra- 
dient in VNP will be useful to supplement additional sci- 
entific data for understanding structure and variation of 
small mammal communities. In this paper, we document 
small mammal community characteristics at VNP based 
on our recent collecting program, and specifically (1) com- 
pare small mammal abundance, diversity and distribution 
in eight habitats (Bamboo, Hagenia Woodland, Brush 
Ridge, Mixed Forest, Herbaceous, Sub Alpine, Alpine 
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Fig. 1. A map of the Virunga massif illustrating the major vegetation zones and the locations of the trapping sites in the VNP. 


zones, and Swamp); (2) assess levels of diversity in rela- 
tion to altitude, environmental variables and habitat char- 
acteristics. (3) Finally, we provide data from two previ- 
ous studies that contribute to assessing changes in Virun- 
ga small mammal diversity and abundance over the past. 


MATERIAL AND METHODS 
The study area 


The Volcanoes National Park, Rwanda (between 
1°21’—1°35’S, 29?22'-29?44' E) is a mountainous region, 
ranging in altitude from 2300 m to 4500 m (Weber 1987). 
It is located on the eastern edge of the Albertine Rift, 
which forms part of the watershed between the Nile and 
Congo River systems (Weber 1987). It lies in north-west- 
em Rwanda and borders its sister parks, the Virunga Na- 
tional Park (ViNP) in the Democratic Republic of the Con- 
go (DRC) and Mgahinga National Park (MGNP) in Ugan- 
da. The Volcanoes National Park (VNP), VINP and MGNP 
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(Fig. 1) form the Virunga volcanoes area (Virunga mas- 
sif). This study was carried out at three sites in the Vol- 
canoes National Park, Rwanda (Mts. Visoke or Bisoke, 
Sabyinyo or Sabinio and Gahinga) from October 2 to No- 
vember 8, 2009. The selection of study sites was based 
on altitude and the major vegetation zones of VNP (Fig. 
1). The major vegetation zones surveyed include: (1) 
Mixed Forest, (2) Bamboo, (3) Hagenia Woodland, (4) 
Herbaceous, (5) Brush Ridge, (6) Sub Alpine zone, and 
(7) Alpine (Plumptre 1991); swamps were categorized as 
an 8th zone. 

For Table 1, we add data from two historical collections. 
The first reflects the efforts of the Swedish Zoological Ex- 
pedition to the ‘Birunga Volcanoes’ from February to 
March, 1921 (Gyldenstolpe 1928). The second collection 
includes the species records from a more recent report 
(Kerbis Peterhans & Austin 1996) that gathered baseline 
data on terrestrial small mammal communities surround- 
ing Kabiranyuma Swamp on the saddle between Muhavu- 
ra and Mgahinga (2980 m) as well as land reclaimed for 
the park (2810 m) from evicted settlers, just on the Ugan- 
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Table 1. Small mammal species recorded in the Virunga massif then and now. Area codes: Burunga (B), Chahafi (Ch), Karisim- 
bi (Ka), Kibati (Ki), Lulenga (L), Lake Mutanda (LM), Mgahinga (Mg), Muhavura (Mu), Mikeno (Mi), Ninagongo (Ni), Ngoma 
(Ng), Sabinio(Sa), Tamohanga (Ta), Tsitsilonga (Ts), Visoke (Vi). References (Ref.): (1) Gyldenstolpe (1928); (2) Kerbis Peter- 


hans & Austin (1996); (3) this study. 
Species 


Crocidura niobe Thomas, 1906 
Crocidura olivieri kivu Osgood, 1910 


Crocidura tarella Dollman, 1915 


Paracrocidura maxima Heim de Balsac, 1959 
Suncus megalura (Jentink, 1888) 
Sylvisorex lunaris ruandae Lonnberg & Gyldenstolpe, 1925 


Sylvisorex vulcanorum Hutterer & Verheyen, 1985 
Myosorex babaulti Heim de Balsac & Lamotte, 1956 
Heterohyrax helgei Lónnberg & Gyldenstolpe, 1925 
Hystrix stegmanni Muller, 1910 

Aethosciurus ruwenzorii vulcanius Thomas, 1909 
Funisciurus carruthersi birungensis Gyldenstolpe, 1927 


Tamsicus vulcanorum vulcanorum Thomas, 1918 
Graphiurus murinus vulcanicus Lonnberg & Gyldenstolpe 1925 


Delanymys brooksi Hayman, 1962 
Dendromus c.f. insignis kivu Thomas, 1916 


Lophuromys aquilus laticeps Thomas & Wroughton, 1907 


Lophuromys woosnami Thomas, 1906 


Arvicanthis abyssinicus rubescens Wroughton, 1909 
Dasymys incomtus medius Thomas, 1906 

Dasymys c.f. rwandae 

Grammomys c.f. dolichurus 

Hylomyscus aeta weileri Lónnberg & Gyldenstolpe, 1925 
Hylomyscus vulcanorum Lónnberg & Gyldenstolpe, 1925 


Lemniscomys striatus cf. massaicus (Pagenstecher, 1885) 
Mus bufo bufo Thomas, 1906 


Mus gratus gratus Thomas & Wroughton, 1910 
Mus triton birungensis Lónnberg & Gyldenstolpe, 1925 


Oenomys hypoxanthus (Pucheran, 1855) 


Praomys degraaffi Van der Straeten & Kerbis Peterhans, 1999 


Praomys jacksoni montis Thomas & Wroughton, 1910 
Thamnomys kempi Dollman, 1911 


Otomys denti kempi Dollman, 1915 
Otomys tropicalis vulcanicus Lónnberg & Gyldenstolpe, 1925 


Tachyoryctes ruandae Lónnberg & Gyldenstolpe, 1925 
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da side of the Virunga Volcanos (Mgahinga Gorilla Na- 
tional Park). Authors of taxon names of small mammals 
are included in Table 1 or may be found in Wilson & Reed- 
er (2005). 


Field methods 


Both live and snap traps were used. The live traps were 
large Sherman folding traps (3"x3"x10"), and small, non- 
folding perforated traps (2”x 2 '2”5 x 6- 2”). Both Live 
and Sherman traps were laid out in lines, 10 m apart (Tews 
et al. 2004), with five lines located randomly in each habi- 
tat at least 100 m away from human paths (Sutherland 
2008). Snap traps were placed at the same points as live 
traps. Each trap was set for three nights. All traps were 
baited with crushed ground-nuts, meat and fish, bananas, 
and pieces of potatoes. They were checked daily between 
0700h and 1000h and again between 1400h and 1600h. 
A locally made trap was used to collect mole rats (Zachy- 
oryctes ruandae) in private cultivated plots adjacent to the 
park. 

At each station, trapped animals were collected, 
processed and/or released (Nicolas & Colyn 2006). Live 
animals were placed in a zip-lock plastic bag and 
weighed using either 100 g or 500 g Pesola spring bal- 
ances. A field identification was provided, sex determined 
and various measurements taken (head-body length, tail 
length, ear length, and hind foot length) using a ruler grad- 
ed in mm. Each newly captured animal to be released was 
marked by fur-clipping to ensure it could be recognized 
if subsequently recaptured. 

Within each of the habitat types, the major plant species 
were recorded from the grids as follows: herbaceous zones 
used 2m’ plots; shrub areas used 5m’ plots; areas domi- 
nated by trees used 10m’ plots. Plant species were iden- 
tified from herbarium material at Karisoke Research Cen- 
ter (KRC) and from the Flora of Rwanda (Troupin 
1977—1988). Temperature and wind speed were collected 
daily between 9:00-9:30 am, allowing comparison be- 
tween different sites. 


Data analysis 


To standardize data for all habitats sampled so that species 
richness and diversity could be compared, trap success for 
each trapline (the number of animals caught per 100 trap 
nights) was calculated. The number of species trapped (per 
100 trap nights) was also calculated for each site using the 
expression: [Species richness/Trap effort] x 100. The 
Shannon- Wiener (H) index of diversity was calculated for 
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all sites sampled and evenness values were derived from 


It. 
H= J -iPitineiy 


with H = the Shannon diversity index; P; = fraction of the 
entire population made up of species I; S = numbers of 
species encountered; and $, = sum from species 1 to 
species S. The Shannon index increases with the number 
of species in a community, and in theory, the Shannon in- 
dex cannot exceed 5.0 (Krebs 1989). A high value of Shan- 
non-Wiener index may indicate a large number for the 
species with similar abundances; a low value indicates 
lower species richness or domination by a few species. 
Evenness values (E)[E = H / In (S)] indicate how num- 
bers of individuals are distributed among species in a com- 
munity. When the evenness is high, the mammal fauna is 
more diverse and the species are equally abundant (Magur- 
ran 2005). 

We tested non-parametric estimators available in Esti- 
mateS software (available at http://viceroy.eeb.uconn.edu 
/estimates/EstimateS) to estimate species richness as an 
alternative to the observed number of species in model 
Chaol. The classic richness estimators Chaol was com- 
puted along with log-linear 95 % confidence intervals 
(Chao 2005, Colwell & Coddington 1994). 

In order to estimate the number of additional species 
necessary to reach an asymptotic plateau, species accu- 
mulation rates were examined and were compared among 
each habitat type (Colwell et al. 2004). The calculation of 
species accumulation rates and the estimation of asymp- 
totes/plateaux (species richness) were derived from the 
distribution of all individuals caught across habitats; 
species richness was estimated as a function of number 
of samples. If the curves became flatter with increasing 
effort (Chao 1 confidence bounds met or connected), it be- 
comes less likely to detect new species in the additional 
samples. 

Similarities of the small mammal communities in the 
different habitats were assessed using Bray-Curtis simi- 
larity indices and linked Cluster Analysis. Tukey's Hon- 
estly-Significant-Difference (Tukey HSD) tests were 
used for post-hoc comparisons when significant among- 
groups differences were indicated. To investigate the re- 
lationship between small mammal species and environ- 
mental factors, a Canonical Correspondence Analysis 
(CCA) was used. This is an indirect gradient analysis tech- 
nique used to study the distribution of species along eas- 
ily measured, recognizable environmental variables 
(Gauch & Whittaker 1972). The mole rat Tachyoryctes 
ruandae was found in cultivated fields outside the park, 
and it was not included in statistical analyses. 
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Table 2. Total numbers of individuals per species captured in the eight habitat types at VNP. 
Species /Habitat/Altitude Mixed Swamp Bamboo  Hagenia Herbaceous — Brush Sub Alpine Alpine Total 
(m) Forest Woodland Ridge 

2380—2580 2540 2540-2660 2740-3020 2900-3183 3220-3400 3420-3600 3640-3710 
Crocidura tarella ] 1 
Crocidura olivieri 3 1 2 1 2 9 
Galerella sanguinea 1 ] 
Graphiurus murinus 2 1 3 6 
Hylomyscus vulcanorum 1 1 2 3 7 2 3 19 
Lophuromys aquilus 8 15 4 8 24 15 15 10 99 
Lophuromys woosnami 10 7 1] 11 13 21 5 78 
Mus bufo 1 l 2 
Oenomys hypoxanthus l l 
Praomys degraaffi 1 1 
Sylvisorex vulcanorum ] ] 
Tachyoryctes ruandae (2) 
TOTALS 24 24 18 29 47 41 23 12 218 
Trap nights 300 300 300 300 300 300 300 300 2400 
Density (Trap success %) 8.0 8.0 6.0 9.7 15.7 13.7 741 4.0 9.1 
# species 5 4 4 8 5 6 3 2 12 
RESULTS 


Species richness 


Small mammal trapping resulted in the capture of 220 in- 
dividuals (208 rodents, 11 shrews, and one mongoose). 
Ninety-nine individuals were female while 121 were male. 
A total of eight rodent species, three shrews and one mon- 
goose species were recorded. All rodents caught belonged 
to the subfamilies Murinae, Deomyinae, Rhizomyinae and 
Graphiurinae (Wilson & Reeder 2005). The most common 
species were Lophuromys aquilus (n=99) and Lophuromys 
woosnami (n—78), combining for 80.54 % of the individ- 
uals captured. These species were encountered in all habi- 
tats; all other species were rare («19.46 % of individuals 
captured in all habitats). 

Small mammal species diversity was highest in Hage- 
nia Woodland (8), intermediate in Brush Ridge (6), Herba- 
ceous and Mixed Forest (5), Swamp and Bamboo (4) and 
lowest in Sub Alpine (3) and Alpine (2). Rodent captures 
were always more numerous than shrew captures (Table 
2). 

For Sub Alpine and Alpine zone Chaol species richness 
estimator upper bounds reached horizontal plateaux 
(asymptotes), and its upper bound reached the lower bound 
for the bamboo and mixed forest (Fig. 2), so sampling ef- 
fort can be considered adequate to reflect actual species 
richness. In other habitats, plateaux were not attained. Fur- 
thermore, species accumulation rates for habitat types and 
altitudinal range showed that the most species rich habi- 
tats in VNP were mid elevation habitats: Hagenia Wood- 
land followed by Brush Ridge and Herbaceous. 
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Species diversity between habitats 


A Bray-Curtis dendrogram based on all individuals cap- 
tured showed that two species captured in Alpine repre- 
sented distinct communities than those from other habi- 
tats. The Alpine (highest elevation at over 3600 m) and 
Bamboo habitats (lowest elevations of 2540—2660 m) had 
the least similar mammal communities (26.66 95) and 
species richness was low in both zones. 

The effective number of species (Shannon diversity) was 
highest (H>1.0) in Hagenia Woodlands and Brush Ridge, 
intermediate (0.8>H<0.6) in Herbaceous, Mixed Forest, 
Bamboo and Swamp and lowest in the Sub Alpine (H=0.2) 
and Alpine habitats (H=0.1). These differences in diver- 
sity indices were significant (Kruskal-Wallis test: 
H- 49=20.49; p-0.005). The variation in community 
species (Evenness) also differed significantly among habi- 
tats types (Kruskal-Wallis test: H} 4,720.49; p=0.005). The 
variation in community species were generally low (<0.5), 
indicating that communities in all habitats were dominat- 
ed by two species (Lophuromys aquilus and Lophuromys 
woosnami). Species diversity was negatively correlated 
with altitude. Altitude accounted for 20 % of the varia- 
tion in species diversity with increasing altitude (p=0.004). 


Environmental correlates of community structures 


The CCA-biplot (the first two axes) for environmental 
variables, habitat and small mammal abundances ex- 


©ZFMK 


182 Deogratias Tuyisingize et al. 


Mixed forest Swamp 


12 


10 


5 10 15 20 5 10 15 20 
Bamboo Hagenia 


"797974999 


coxa 
eee 


5 10 15 
Herbaceous 


Species richness 
10 12 14 


6 8 


4 


20 -30 


Subalpine 


19 20 


Number of individuals sampled 
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circle represented Chaol estimate, while dotted line is 95% confidence interval. 
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plained 19.3+39.4 (equals to 58.7 %) of variance in the 
data set (Fig. 3). The environmental variables that affect- 
ed separation of sites and species along axis | most strong- 
ly were temperature, altitude, canopy cover, herbaceous 
cover and wind speed. 

Lophuromys aquilus, Lophuromys woosnami and Hy- 
lomyscus vulcanorum plotted near the intersection of the 
axes, suggesting that these species are not strongly influ- 
enced by any of the six environmental variables and that 
they are habitat generalists. Furthermore, Crocidura 
olivieri kivu also seemed to be a generalist, and occurred 
in five habitat types. The percentage of canopy cover has 
a greater influence on the distribution of Graphiurus mur- 
inus vulcanicus, a primarily arboreal species. 


DISCUSSION 


Our study found few small mammal species at VNP com- 
pared to previous studies. The contiguous forest at 
Mgahinga National Park has 16 known rodent species and 
seven shrew species (Kerbis Peterhans & Austin 1996), 
eight of which were recorded during this study: Hy- 
lomyscus vulcanorum; Lophuromys aquilus; Lophuromys 
woosnami; Oenomys hypoxanthus; Sylvisorex vulcano- 
rum; Crocidura olivieri; Graphiurus murinus and Mus bu- 
fo. Differences in size of the protected areas, trapping 
techniques, time of year, elevation, and habitats sampled 
may account for differences in species diversity with this 
adjacent site. The differences between sizes of sampled 
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areas and trapping methods might have influenced differ- 
ences. For instance squirrels can be obtained using shot- 
guns. Shrews are more easily captured with pitfalls. 
Gyldenstolpe (1928) covered the whole Volcanoes Nation- 
al Park, while this study has focussed on Mt Bisoke and 
the foot of Mt Sabyinyo. 

Our study confirms a pattern of variation in species rich- 
ness of small mammal communities along an altitudinal 
gradient. Habitat heterogeneity and altitudinal variations 
are a major factor affecting small mammal diversity (Is- 
abirye-Basuta & Kasenene 1987; Stanley & Hutterer 
2007). Mid elevation peaks in species diversity are com- 
mon in both rodents and shrews whereas increasing ele- 
vation leads to fewer species. Lophuromys aquilus and Lo- 
phuromys wosnami are adaptable to the VNP habitats be- 
cause they require moist scrub, forest, grassy, and tropi- 
cal moist montane areas (Kingdon 1984). The elevation- 
al diversity might be related to climatic impacted factors 
such as rainfall, temperature, productivity, competition, re- 
source abundance, habitat complexity, or habitat diversi- 
ty (Lomolino 2001). 

Species accumulation curves show insufficient sampling 
effort in six of eight studied habitats. Further, extreme 
weather during our sampling periods may have reduced 
trapping success. Some species previously collected in 
VNP were not captured during this study including Cro- 
cidura niobe, Suncus megalura, Sylvisorex lunaris ruan- 
dae, Paracrocidura maxima, Myosorex babaulti and 
eleven mice/rat species including Delanymys brooksi, 
Dendromus insignis kivu, Dasymys incomptus, Dasymys 
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Fig. 3. 


CCA-biplot showing the influence of various environmental variables on small mammal distributions during this study 


period. Gram: Graphiurus murinus; Crot; Crocidura tarella; Croc: Crocidura olivieri; Hylv: Hylomyscus vulcanorum; Lapq: Lo- 
phuromys aquilus; Lopw: Lophuromys woosnami; Musb: Mus bufo; Oenh: Oenomys hypoxanthus; Prad: Praomys degraaffi; Sylv: 


Sylvisorex vulcanorum; Gals: Galerella sanguinea. 
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c.f. rwandae, Grammomys c.f. dolichurus, Hylomyscus ae- 
ta, Thamnomys kempi, Otomys denti kempi, Otomys trop- 
icalis vulcanus, Mus triton, Mus gratus (Gyldenstolpe 
1928; Hutterer et al. 1987; Kerbis Peterhans & Austin 
1996). More trapping over different seasons and with the 
added use of pitfall traps in more habitats are needed to 
adequately document small mammal communities at V NP. 

The study area contains relatively few endemic species 
(six) compared with other Albertine Rift sites. This com- 
pares with twenty-one endemic species in the Ruwenzori 
Mountains of Uganda (Kerbis Peterhans et al. 1998), and 
18 endemic species in Kahuzi-Biega forest (Kaleme et al. 
2007). Five of the six endemic species and the one threat- 
ened species were found at low and middle altitudes at 
VNP. 

Species from VNP have a high probability of survival 
owing to international and governmental protection efforts 
directed at the mountain gorilla (Gorilla beringei 
beringei), an ‘umbrella species’ from a conservation per- 
spective. Since gorillas and endemic rodents and shrews 
depend on the same habitats, the conservation of the first 
acts to ensure the conservation of the second. In order to 
maximize the maintenance of biodiversity conservation in 
general and small mammals in particular, conservation ef- 
forts need to target both low elevation and mid elevation 
zones at VNP. 

Our results have implications for how small mammal 
communities in VNP may alter with climate change. With 
increasing temperature, formerly low-elevation small 
mammal species may expand their ranges upwards while 
those of high-elevation species may contract, leading to 
changed community composition at mid- and high eleva- 
tions (Moritz et al. 2008). More long-term monitoring is 
desirable to incorporate these factors, and to examine in 
detail how to mitigate these risks. 

Our study indicates that the small mammal species have 
changed much since the pioneer study of Gyldenstolpe 
(1928). This may be a signal of concern for the conser- 
vation of this forest, and stresses the importance of reg- 
ular small mammal surveys with close attention to asso- 
ciated habitats. 

The problem of spatial scales is also important; 40 trap 
lines located at three study areas may not have been suf- 
ficient in a protected area of 150 km’ and 40 km of width. 
Long-term studies are needed to better understand the ef- 
fects of climatic events relative to those of local differ- 
ences related to habitat structure. The findings of species 
accumulation rates and comparatively low richness dur- 
ing this study indicate that the small mammal species list 
for VNP is still incomplete and further inventory work is 
needed. We recommend that future, longer-term studies 
are carried out to confirm our species list, and to exam- 
ine demographic processes and the impact of climatic 
changes on gorilla habitat. 
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